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ABSTRACT 

Combining galaxy cluster data from the ROSAT All-Sky Survey and the Chandra X- 
ray Observatory, cosmic microwave background data from the Wilkinson Microwave 
Anisotropy Probe, and galaxy clustering data from the WigglcZ Dark Energy Survey, 
the 6-degree Field Galaxy Survey and the Sloan Digital Sky Survey III, we test for 
consistency the cosmic growth of structure predicted by General Relativity (GR) and 
the cosmic expansion history predicted by the cosmological constant plus cold dark 
matter paradigm (ACDM). The combination of these three independent, well stud- 
ied measurements of the evolution of the mean energy density and its fluctuations 
is able to break strong degeneracies between model parameters. We model the key 
properties of cosmic growth with the normalization of the matter power spectrum, 
o"8i an d the cosmic growth index, 7, and those of cosmic expansion with the mean 
matter density, fi m , the Hubble constant, Hq, and a kinematical parameter equiva- 
lent to that for the dark energy equation of state, w. For a spatially flat geometry, 
w = —1, and allowing for systematic uncertainties, we obtain cr 8 = 0.785 ± 0.019 and 
7 = 0.570^0 063 ( a t the 68.3 per cent confidence level). Allowing both w and 7 to vary 
we find w = -0.950±g;{$g and 7 = 0.533 ± 0.080. To further tighten the constraints 
on the expansion parameters, we also include supernova, Cepheid variable and baryon 
acoustic oscillation data. For w = —1, we have 7 = 0.616±0.061. For our most general 
model with a free w, we measure tt m = 0.278t°;^i, H Q = 70.0 ± 1.3 km s^Mpc -1 

and w = — 0.987lo'o53 for the expansion parameters, and er 8 = 0.789 ± 0.019 and 
7 = 0.604 ±0.078 for the growth parameters. These results are in excellent agreement 
with GR+ACDM (7 ~ 0.55; w = —1) and represent the tightest and most robust 
simultaneous constraint on cosmic growth and expansion to date. 

Key words: cosmological parameters - cosmology: observations - X-ray: galaxies: 
clusters - cosmic microwave background - large-scale structure of Universe 



1 INTRODUCTION 

The unexpected measurement from type la supernova 
(SNIa) data of late-time cosmic acceleration by Riess et al. 
(1998) and Perlmutter et al. (1999) initiated a series of 
theoretical and observational efforts to unveil the na- 
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ture of its underlying cause. However, to this day it is 
still unclear whether the origin of this phenomenon is 
due to a new energy component, spurious cosmological 
assumptions, or modifications of gravity at large scales. 
A number of theoretical approaches and observational 
probes have been developed to investigate these different 
possibilities (for recent reviews see Copeland et al. 2006; 
Frieman et al. 2008; Allen et al. 2011; Clifton et al. 2012; 
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Weinberg et al. 2012). Current data on the energy content, 
geometry, and expansion and growth histories of the Uni- 
verse do not show any deviation from the standard cos- 
mological paradigm, ACDM (Allen et al. 2008; Mantz et al. 
2008, 2010a; Vikhlinin et al. 2009; Percival et al. 2010; 
Reid et al. 2012; Conley et al. 2011; Blake et al. 2011a,b, 
2012; Suzuki et al. 2012; Hinshaw et al. 2012). However, 
the cosmological constant model suffers from well-known, 
serious theoretical problems that present-day dark en- 
ergy models have not been able to improve upon. 
For modified gravity models, various approaches have 
been developed 1 : parameterized frameworks (Hu & Sawicki 
2007; Bertschinger & Zukin 2008; Amin et al. 2008), con- 
sistency tests of GR (Linder 2005; Linder & Cahn 
2007; Di Porto & Amendola 2008; Zhang et al. 2007; 
Nesseris & Perivolaropoulos 2008; Acquaviva et al. 2008), 
and alternative theories of gravity (Dvali et al. 2000; 
Carroll et al. 2004; Arkani-Hamed et al. 2004; Nicolis et al. 
2009; de Rham et al. 2011). Recent works have used 
a variety of experiments and data sets to constrain 
gravity properties and models and found no signifi- 
cant deviations from GR (see e.g. Schmidt et al. 2009; 
Rapetti et al. 2009, 2010; Reyes et al. 2010; Daniel et al. 
2010; Zhao et al. 2010, 2012; Giannantonio et al. 2010; 
Wojtak et al. 2011; Hojjati et al. 2011; Lombriser et al. 
2012; Hudson & Turnbull 2012; Basilakos & Pouri 2012; 
Samushia et al. 2013; Simpson et al. 2013). To further test 
the overall standard paradigm, GR+ACDM, it is crucial to 
use data sets able to robustly constrain the key properties 
of the model, and to combine complementary data sets to 
break the degeneracies between the model parameters. 

Rapetti et al. (2010, hereafter R10) tested GR using 
ROSAT All-Sky Survey (RASS) and Chandra X-ray Ob- 
servatory (CXO) data of cluster abundance and scaling re- 
lations (Mantz et al. 2010a,b, hereafter M10a,b). R10 ob- 
tained strong constraints on GR, even when marginaliz- 
ing over conservative systematic and astrophysical mod- 
eling uncertainties in the evolution of the cluster X- 
ray luminosity-mass relation. When combining the clus- 
ter growth data with measurements of the anisotropies in 
the cosmic microwave background (CMB) from the Wilkin- 
son Microwave Anisotropy Probe ( WMAP; Spergel et al. 
2003, 2007; Komatsu et al. 2009; Dunkley et al. 2009; 
Komatsu et al. 2011; Hinshaw et al. 2012; Bennett et al. 
2012, and companion papers), they highlighted a large de- 
generacy between 7 and ag, which limited the constraints on 
each parameter individually. Here we include complemen- 
tary data that break this degeneracy. In particular, we use 
measurements of the growth rate and the Hubble parame- 
ter from joint redshift space distortions (RSD) and Alcock- 
Paczynski (AP) effect constraints (Blake et al. 2011c, here- 
after Bll) from the WiggleZ Dark Energy Survey (Wig- 
gleZ; Drinkwater et al. 2010). We also use a low redshift 
RSD constraint (Beutler et al. 2012) from the final data re- 
lease (DR3) of the 6-degree Field Galaxy Survey (6dFGS; 
Jones et al. 2009) and an RSD and AP constraint from the 
latest data release (DR9) of the Sloan Digital Sky Survey 



1 Similar approaches can also be used to study clustering dark 
energy models. 



(SDSS) III Baryon Oscillation Spectroscopic Survey (BOSS; 
Reid et al. 2012). 

In addition to our primary data sets, and to tighten the 
constraints on the expansion parameters, we also present 
results including data from the Union II SNIa sample 
(Suzuki et al. 2012), baryon acoustic oscillation (BAO) 
measurements from a combined analysis (Percival et al. 
2010) of 2-degree Field Galaxy Redshift Survey (2dFGRS; 
Colless et al. 2003) and SDSS-II DR7 (Abazajian et al. 
2009) data as well as from a recent analysis (Reid et al. 
2012) of SDSS-III BOSS data at a higher redshift, and H 
measurements from the Supernovae and Ho for the Equation 
of State program (SHOES; Riess et al. 2011). 

We find that by combining cluster, CMB and galaxy 
data we are able to break the key degeneracies between 7 
and erg and obtain tight and robust constraints on cosmic ex- 
pansion and growth. We model the expansion primarily with 
f2 m , Ho, and w and the growth with cr& and 7. We find that, 
individually, the CMB and galaxy data have large degen- 
eracies in the growth plane but that, crucially, these degen- 
eracies are nearly orthogonal. The individual and combined 
constraints from cluster, CMB and galaxy data are consis- 
tent with one another, making this a very robust measure- 
ment, and in good agreement with GR and ACDM. While 
individually clusters provide the tightest constraints in the 
growth plane, the combination of clusters, the CMB and 
galaxies provides significantly improved constraints and ar- 
guably the most robust measurement of cosmic structure 
growth to date. 



2 COSMOLOGICAL MODEL 

We adopt a purely phenomenological model to conveniently 
test the consistency of current observations with both the 
cosmic expansion history and the cosmic growth history pre- 
dicted by ACDM+GR. 

Our model assumes neither the existence of a new com- 
ponent, dark energy, nor a modification of GR. Instead, the 
parameters of the model represent departures from key kine- 
matical and dynamical features of ACDM+GR. Deviations 
from such benchmarks would indicate disagreement of the 
observed evolution of the background and density perturba- 
tions with the standard cosmological paradigm. 2 



2.1 Cosmic expansion history 

We model the expansion history using the evolution param- 
eter E(a) = H(a)/Ho, where H{a) is the Hubble parameter 
as a function of the scale factor a and Ho its present-day 
value. We parameterize E(a) as follows 

E(a) = [f2 m a' 3 + (1 - fi m ) a- 3(1+ll,) ] V * . (1) 



2 Note that different physical scenarios can cause similar de- 
partures from this paradigm. For example, specific models of 
clustering or interacting dark energy and of modified gravity 
might provide similar deviations from the density perturbations 
of ACDM+GR. 
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f2 m is the present, mean matter density in units of the crit- 
ical density of the Universe and w a kinematical parameter 
inspired by the dark energy equation of state. Since for our 
test we do not assume any particular scenario for cosmic 
acceleration, such as dark energy, we use w only to con- 
veniently fit expansion history data, matching ACDM for 
w = — 1. Below, we present results for two expansion mod- 
els, w = -1 (ACDM) and w constant (tuCDM). For both 
cases, we assume a spatially flat geometry (i.e., the curva- 
ture energy density fi^ = 0) 3 . 



2.2 Cosmic growth and cluster abundance 

We model the growth history at late times by parameteriz- 
ing the linear growth rate of density perturbations on large 
scales, /(a), as a power law of the evolving mean matter 
density, fi m (a) = fi m a~ 3 -E(a) -2 , such as (Peebles 1980; 
Wang & Steinhardt 1998; Linder 2005) 



f(a) 



d\n6 
din a 



n m (a) 7 , 



C2) 



where 7 is the growth index 4 , for which we recover GR when 
7 ~ 0.55. J 8 = 8pm/ pm is the ratio of the comoving matter 
density fluctuations, 8p m , with respect to the cosmic mean, 
pm- While at early times we assume GR, for z < z t we 
obtain 5{z) from equation 2 using as an initial condition 
8(zt) calculated within GR. Normalizing S(z) to S(z t ), we 
obtain the growth factor, D(z) = S(z)/S(z t ). Here we use 
z t = 30, which is well within the dark matter dominated 
era, when f(a) ~ 1 for both the 7-model (equation 2) and 
GR. We then calculate the matter power spectrum of such 
fluctuations for a given wavenumber, k, as 



where Wm (k) is the Fourier transform of the window func- 
tion. From this expression, erf is defined as the z = vari- 
ance in the density field at scales of 8/i _1 Mpc, where as 
is widely used as a parameter for the normalization of the 
matter power spectrum. 

Here we use a(M, z) to calculate the abundance of dark 
matter halos as a function of mass and redshift 



n(M,z) = J 



d In a 
dM' 



■ dM' 



(5) 



where J-(cr, z) is a convenient fitting formula obtained 
from large N-body simulations of dark matter particles 
(Tinker et al. 2008), 



T{a, z) = A 



1 



(6) 



The parameters of this formula have a generic redshift de- 
pendence of the form x(z) — xo(l + z) ea * , with x represent- 
ing A, a, b or c. The values for each xo and « x are given in 
Tinker et al. (2008). As in MlOa, we introduce an additional 
parameter, e, to account for residual systematic uncertain- 
ties in the evolution of J- (a, z) due to non-ACDM scenar- 
ios. Remarkably, J- (a, z) encapsulates the non-linear cosmic 
growth history and appears to be almost universal for a wide 
range of cosmologies (see R10 for more details). 

We marginalize over the uncertainties in the parameters 
of J- (a, z), accounting for their covariance and for additional 
systematic uncertainties due to e.g. the presence of baryons 
following the method described in MlOa. Note, though, that 
as shown in MlOa the uncertainties in J- (a, z) are subdom- 
inant in the analysis. R10 also verified that e is essentially 
uncorrelated with 7. 



P(k,z) oc k ns T 2 (k,z t )D(z) 2 , (3) 

where T(k, z%) is the matter transfer function of GR in the 
synchronous gauge at redshift zt and n s the primordial scalar 
spectral index. 

The variance of the linearly evolved density field, 
smoothed by a spherical top-hat window function of comov- 
ing radius R enclosing mass M = 4-7rp m J? 3 /3, is 

a 2 (M, z) = — J k 2 P(k, z)\W M {k)\ 2 dk , (4) 



3 Using cluster, CMB and SNIa data, Rapetti et al. (2009) found 
a negligible correlation between f2 k and 7. They also showed that 
the constraints on 7 were not significantly weaker when including 

as a free parameter. Note also that if were included as a free 
parameter, an extension of equation 2 proposed by Gong et al. 
(2009) would fit better the predictions from GR. 

4 Many models of modified gravity predict a growth index that 
varies with time and length scale, 7(0, k). Note again, though, 
that here we do not use this parameter as a diagnostic of the true 
theory of gravity, but rather as a consistency test for GR. 

5 For current results, this value is a good approximation to be 
used as a GR reference. At higher accuracy, though, the growth 
index of GR has relatively small redshift and background param- 
eter dependencies (see e.g. Polarski & Gannouji 2008). 



2.3 Integrated Sachs- Wolfe effect 

In our CMB analysis we include the constraint on 7 from the 
Integrated Sachs- Wolfe (ISW) effect of the CMB using the 
method and assumptions described by Rapetti et al. (2009, 
2010). In brief, the low multipoles of the CMB are sensi- 
tive to the growth of cosmic structure due to the effect of 
the time- varying gravitational potentials of large scale struc- 
tures on the CMB photons crossing them. We calculate the 
contribution of these photons to the temperature anisotropy 
power spectrum as (Weller & Lewis 2003) 

Aj sw (fc) = 2 J dte~ T{t U'ji [k(t - t )] , (7) 

where t is the conformal time and to its present-day value, 
r the optical depth to reionization, ji(x) the spherical 
Bessel function for the multipole I, and <j>' the confor- 
mal time variation of the gravitational potential. Taking 
the derivative of the Poisson equation with respect to t, 
we calculate the latter quantity for the 7-model 6 as <fi' — 
AixGa 2 k~ 2 H 8p m [1 — m (a) 7 ] , where H is the conformal 
Hubble parameter. Since the ISW effect is only relevant for 

6 Since here we are testing GR, we assume no contributions to 
<j> from the anisotropic stress and energy flux of the Weyl tensor 
(Challinor & Lasenby 1999). 
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z < 2, as an initial condition to solve this equation we match 
A\ sw (k) to that of GR at z t = 2. 7 

Note, however, that the constraining power on 7 from 
the ISW effect is small compared to that of the cluster data 
(Rapetti et al. 2009). For the current analysis, the primary 
relevance of the CMB is its ability to tightly constrain the 
combination of growth parameters erg and 7 (see Section 5). 

2.4 The Alcock-Paczynski effect and 
redshift-space distortions 

The Alcock-Paczynski test is a geometrical means of prob- 
ing the cosmological model by a comparison of the ob- 
served tangential and radial dimensions of objects which 
are assumed to be isotropic in the correct choice of model. 
It can be applied to the 2-point statistics of galaxy clus- 
tering if the redshift space distortions, the principal addi- 
tional source of anisotropy, can be successfully modelled 
(Ballinger et al. 1996; Matsubara & Suto 1996; Matsubara 
2000; Seo & Eisenstein 2003; Simpson & Peacock 2010). By 
equating radial and tangential physical scales, the AP test 
determines the observable F(z) = (1 + z)Da(z)H(z)/c, 
where Da(z) is the physical angular diameter distance and 
c is the speed of light. 

In the model fit for F(z), the normalized growth rate, 
fag(z), is determined simultaneously. Here f(z) is again the 
logarithmic rate of change of the growth factor at redshift 
z (see equation 2) and cr g (z) = [D(z)/D(0)] <7g. In Bll, 
RSD were modelled using the fitting formulae provided by 
Jennings et al. (2011) to determine the density-velocity and 
velocity-velocity power spectra, marginalizing over a linear 
bias factor. Tests were performed to ensure that the re- 
sults were not very sensitive to the model used for the non- 
linear RSD, the real-space power spectrum, or the range 
of scales fitted (fc max < 0.2 /iMpc -1 , for the measurements 
used here). 

For a low-redshift survey such as 6dFGS, the Alcock- 
Paczynski distortion is negligible (since distances in 
h^ 1 Mpc are approximately independent of the assumed cos- 
mological model) 8 . For 6dFGS, the growth rate measure- 
ment of Beutler et al. (2012) was obtained by again assum- 
ing the model of Jennings et al. (2011) to described non- 
linear RSD. 

For the BOSS measurements of the RSD and AP effect, 
the modeling of the matter density and velocity fields was 
performed following the approach of Reid & White (2011). 
The latter uses perturbation theory to calculate the non- 
linear redshift space clustering of halos in the quasilinear 
regime and the halo model framework to describe the galaxy- 
halo relation. This model was tested against a large set of 
galaxy catalogs from N-body simulations and only fit over 
those scales where the quasilinear velocity field was thought 
to dominate the signal and the small-scale random velocities 
could be simply modeled and marginalized over. 

For all the RSD and AP effect measurements employed 

7 For our analysis, the difference from calculating 8(z) using zt 
equals to 2 or 30 is negligible since both redshifts are well within 
the dark matter dominated era, when f(a) tends to 1 for any 7. 

8 Beutler et al. (2012) calculated the uncertainties in F(z) and 
showed that they are unimportant. 



in the paper (see Section 4.2), the parameters used to fit the 
2D galaxy power spectrum and galaxy correlation function 
data have negligible covariance with the parameters in the 
current analysis. Also, the linear model as well as the non- 
linear corrections assumed in those analyses lie within the 
GR+ACDM paradigm tested here 9 . 



3 PHYSICS OF THE OBSERVABLES 

In this section, we describe the physical mechanisms behind 
the principle degeneracies between our most relevant growth 
and expansion parameters, for each of our primary observa- 
tions. 

3.1 CMB anisotropics 

From the normalization and tilt of the CMB temperature 
anisotropy power spectrum, we can primarily constrain the 
scalar amplitude and spectral index of primordial fluctua- 
tions; from the position of its first peak, the mean energy 
density of curvature and dark energy; and from the ampli- 
tudes of the second and third peaks, those of dark and bary- 
onic matter. These measurements provide strong constraints 
on the content of the background energy density and its lin- 
ear density fluctuations at high redshift. For a given value 
of the growth index, 7, these translate into tight constraints 
on the amplitude of the matter power spectrum today, as- A 
model with faster perturbation growth, i.e. with a small 7, 
implies large fluctuations today, i.e. large erg, and vice- versa. 
This provides a large, negative correlation between erg and 
7 (see Figure 1). At low redshift, the ISW effect of the CMB 
data (see Section 2.3) constrains 7, which is otherwise un- 
constrained by this data set. 

3.2 Distribution of galaxies 

From measurements of the anisotropic clustering of galax- 
ies, we use constraints on the product f(z) erg (z) and on the 
quantity F(z), where the latter are purely expansion his- 
tory constraints, i.e. on f2 m (z). Both of these constraints, 
from RSD and AP effect measurements respectively, are re- 
quired to measure 7 = In f(z)/ In f2 m (z) from galaxy data 
alone. 10 The current uncertainty on the linear galaxy bias, 
b(z), limits the ability to measure as from the normalization 
of the galaxy power spectrum, which scales with ag(z) b(z), 
and to measure f(z) using RSD constraints on f(z)/b(z), as 
previously commonly used. As proposed by Song & Percival 
(2009), here we use instead RSD constraints on f(z)as(z), 
which are independent of b(z), and obtain a positive corre- 
lation between 7 and as (see Figure 1) for a ACDM ex- 
pansion model and data within a relatively low-z range, 
where f(z) increases towards 1. The faster the perturba- 
tions grow (small 7), the smaller the present-day perturba- 
tion amplitude, erg, needs to be to provide the same amount 

9 The non-linear modeling from Jennings et al. (2011) used in 
the WiggleZ and 6dFGS analyses also encompasses a range of 
quintessence dark energy models. 

10 In the same way as for the AP effect, the addition of the BAO 
constraints on f2 m (,z) improves significantly the measurement of 
7 for the combination gal+BAO (see the right panel of Figure 1). 
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of anisotropy in the distribution of galaxies at redshift z. At 
higher-z, where /(z) ~ 1 and f(z) as(z) ~ os(z), the corre- 
lation between 7 and as becomes negative (see Section 3.3). 
Adding high-z data from future missions will then help to 
break the large degeneracy of the current data between 7 
and <rg. 

3.3 Cluster abundance and masses 

For clusters, we have direct constraints on as(z) and fi m (z) 
from abundance, mass calibration and gas mass fraction 
data (see Sections 2.2 and 4.1). ag(z) measurements pro- 
vide us with constraints not only on a$(z = 0), from the 
local cluster mass function, but also on the growth rate 
/(z) — —(1 + z)d\nas(z)/dz. Together, the constraints on 
as{z) and fi m (z) constrain 7. 

The evolution of ug(z) = <78e~ 9 ' z ' depends on 7, Q m 
and w as follows 

g(z) = r(l + /)- 1 [p(/)-l]" 7 p(/r^' (8) 
Jo 

= (Swjy 1 [A(z) - A(0)] , (9) 

where A(z) = \p(z) - l] 1_7 p(2) 7 aft [1, 1; 1 + 75P(*)]. a^i 
is a hypergeometric function, p(z) = po(l + z)~' iw and 
po = fim/tfim — 1). In practice, a negative degeneracy be- 
tween as and 7 exists due to the limited precision of clus- 
ter mass estimates, but it is notably smaller than those de- 
scribed above (see Figure 1). Within the precision of the 
data, indistinguishable cluster count evolution can be pro- 
duced by models with e.g. as of 0.8 and a growth rate con- 
sistent with GR, or with a slightly larger present-day ampli- 
tude and faster growth (smaller 7), for which crg(z) decreases 
with z a bit more steeply. 

For the 7+iuCDM model, the dependence of as(z) on 
the product w 7 implies a negative correlation on the w, 7 
plane (see Figure 2). Within the precision of the data, a 
fast expansion history (small w) can be mimicked by a slow 
growth history (large 7), and vice- versa. 



4 DATA ANALYSIS 
4.1 Galaxy cluster data 

For clusters we use two experiments: growth of structure 
(M10a,b) and gas mass fraction (/ gas ; Allen et al. 2008) 1 . 

Following the methods developed by M10a,b for the 
cluster growth analysis, we self-consistently and simultane- 
ously combine X-ray survey and follow-up observations to 
obtain the best constraints possible while accounting fully 
for selection biases. We employ the survey data to deter- 
mine cluster abundances and the follow-up data to calibrate 
cluster masses from three observables: luminosity, tempera- 
ture and X-ray emitting gas mass. For the survey data we 
employ three wide-area cluster samples drawn from RASS: 
the Bright Cluster Sample in the northern sky (BCS; z < 
0.3 and F x (0.1 -2.4keV) > 4.4 x 10" 12 ergs -1 cm" 2 ), the 

11 Note that the cluster growth analysis employs the / gas analysis 
to calibrate the masses for the scaling relations of Section 4.1.1 
using gas mass as a proxy for total mass (sec details in MlOa). 



ROSAT-ESO Flux Limited X-ray sample in the southern sky 
(REFLEX; z < 0.3 and F x > 3.0 x 10~ 12 ergs" 1 cm" 2 ), and 
the Bright Massive Cluster Survey with ~ 55 per cent sky 
coverage (Bright MACS; 0.3 < z < 0.5 and F x > 2 x 10" 12 
ergs" 1 cm" 2 ). To keep systematic uncertainties to a mini- 
mum and maintain a trivial constant scaling between X-ray 
gas mass and total mass, for all three samples we impose a 
lower luminosity cut of 2.5 x 10 44 /if 2 ergs" 1 (0.1 - 2.4 keV) 
leaving a total of 78 clusters from BCS; 126 clusters from 
REFLEX; and 34 clusters from Bright MACS. In total we 
use 238 clusters. For 94 of these clusters, we employ follow- 
up observations from CXO or pointed observations from 
ROSAT (MlOb; distributed along the same redshift range 
of the survey data < z < 0.5) to constrain simultaneously 
the luminosity-mass (L-M) and temperature-mass (T-M) 
relations using the model from MlOb (see a brief description 
in Section 4.1.1). 

For the / ga a analysis, we use the methods and data set 
of Allen et al. (2008) for 42 massive, hot (kT > 5keV), dy- 
namically relaxed, X-ray luminous galaxy clusters spanning 
the redshift range 0.05 < z < 1.1. 

4-1.1 Scaling relations model 

We model the L-M scaling relation as (MlOb) 

(t(m)) = d m + p{ m m + P l 2 m log 10 (l + z) , (10) 
with a log-normal intrinsic scatter at a given mass of 

aim(z) — a em (l + a' (m z) , (11) 

where I = log 10 [L5oo-E(z) _1 /10 44 ergs" 1 ] and m = 
log 10 [£(z)M5oo/10 15 M Q \. The subscript 500 refers to quan- 
tities measured within radius rsoo, at which the mean, en- 
closed density is 500 times the critical density of the Universe 
at redshift z. We model the T-M scaling relation (t(m)), 
where t = log 10 (fcTsoo/ keV), and its scatter at m (z) us- 
ing the same equations 10 and 11 but with the parameters 
/3o"\ Pi™, Pi™, °~tm and a' tm instead of those with index I. 
When /3| m = and = we have "self-similar" evolution 
of the L-M and T-M relations respectively (Kaiser 1986; 
Bryan & Norman 1998) 12 . a' lm — and a' tm = correspond 
to scaling relations with non-evolving scatter. 

MlOb showed that current data do not require depar- 
tures from self-similar evolution and constant scatter. R10 
demonstrated that 7 correlates weakly with departures from 
self-similarity and constant scatter in the L-M relation and 
negligibly for those in the T-M relation. Here we therefore 
assume self-similar evolution and constant scatter for both 
relations (/?| m = a' lm = fe m = a' tm = 0). 

4.2 Galaxy clustering data 

For WiggleZ, a series of growth and expansion analyses have 
recently been released, and here we build on one in particu- 
lar: the joint analysis of the AP effect and growth of struc- 
ture presented by Bll, which contains four redshift bins of 

12 Self-similar evolution is entirely determined by the E(z) fac- 
tors in the definitions of I, t and m. 
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width Az = 0.2, spanning the redshift range 0.1 < z < 0.9. 
The WiggleZ survey at the Australian Astronomical Obser- 
vatory was designed to extend the study of large-scale struc- 
ture over large cosmic volumes to higher redshifts z > 0.5, 
complementing SDSS observations at lower redshifts. The 
survey, which began in August 2006, completed observations 
in January 2011 and has obtained of order 200,000 redshifts 
for UV-bright emission-line galaxies covering of order 1000 
square degrees of equatorial sky. 

For the WiggleZ analysis we fit our cosmological mod- 
els to the joint measurements of RSD and AP distortion 
presented by Bll. For this, we use the constraints ob- 
tained by Bll as a bivariate Gaussian likelihood for / ag(z) 
and F(z), including the large correlations between them. 
From Bll, we have four bins with effective redshifts z = 
(0.22,0.41,0.60,0.78) and f a 8 (z) = (0.53 ± 0.14, 0.40 ± 
0.13,0.37 ± 0.08,0.49 ± 0.12), F(z) = (0.28 ± 0.04, 0.44 ± 
0.07,0.68 ± 0.06,0.97 ± 0.12) and correlation coefficients 
r = (0.83,0.94,0.89,0.84). 

For the 6dFGS analysis we use the growth rate of struc- 
ture measurement obtained by Beutler et al. (2012). The 
6dFGS is a combined redshift and peculiar velocity survey 
covering nearly the entire southern sky with the exception of 
a 10 degree band along the Galactic plane. Observed galax- 
ies were selected from the 2MASS Extended Source Catalog 
(Jarrett et al. 2000) and the redshifts were obtained with the 
6-degree Field multi-fibre instrument at the U.K. Schmidt 
Telescope between 2001 and 2006. The final 6dFGS sample 
contains about 125,000 galaxies in 5 bands distributed over 
~ 17,000 square degrees with a mean redshift of z — 0.052. 

For the analysis of the RSD from 6dFGS data we use the 
constraints obtained by Beutler et al. (2012) as a Gaussian 
likelihood for / ag(z) = 0.423 ±0.055 at an effective redshift 
z = 0.067. 

The analysis of the SDSS-III BOSS results from 
Reid et al. (2012) are based on the high-z sample CMASS, 
which consists of 264,283 galaxies in the redshift range 
0.43 < 2 < 0.7 over 3,275 square degrees. As part of SDSS- 
III (Eisenstein et al. 2011), BOSS has imaged the South 
Galactic sky for an additional 3100 square degrees over 
SDSS-II. This has increased the total sky coverage of SDSS 
imaging to 14,055 square degrees. As its primary goal, BOSS 
targets for spectroscopy luminous galaxies selected from the 
SDSS imaging. Within BOSS, CMASS is a roughly volume- 
limited sample of massive, luminous galaxies (for more detail 
see e.g. Masters et al. 2011) tracing a cosmological volume 
at a high enough density to enable powerful statistical stud- 
ies of large-scale structure. 

For the analysis of the growth rate and AP effect mea- 
surements of CMASS BOSS, we use a bivariate Gaussian 
likelihood for f <r 8 (z) = 0.43 ± 0.07 and F(z) = 0.68 ± 0.04 
with a correlation coefficient r — 0.87 at an effective redshift 
z = 0.57 (Reid et al. 2012) 13 . Note that this redshift is sim- 
ilar to that of the third redshift bin of the WiggleZ analysis, 
z — 0.6. Due to the small overlap and the uncorrelated shot 
noise between the two surveys, their covariance should be 



3 For the results in Section 5 that include these and the distance- 
scale constraints from the BAO signature in the CMASS BOSS 
data, wc extend this likelihood to account for the correlations 
between these three measurements as discussed in Section 4.4. 



minimal. Importantly, the results obtained by the two in- 
dependent experiments, which target very different galaxy 
types, and require very specific studies of their nonlinear 
properties and modeling uncertainties, are consistent. 

4.3 CMB data 

For the CMB analysis, we use the data and likelihood code 14 
from WMAP 15 . For the analyses including CMB data, we 
also fit for the mean physical baryon and dark matter densi- 
ties, flbh 2 and Q c h 2 , the optical depth to reionization, r, the 
logarithm of the adiabatic scalar amplitude, ln(A s ), which 
is related to a$, and the adiabatic scalar spectral index, n s . 
For these analyses, instead of Ho we fit 9, the (approximate) 
ratio of the sound horizon at last scattering to the angular 
diameter distance, which is less correlated with other param- 
eters than Ho (Kosowsky et al. 2002). We also marginalize 
over the amplitude of the Sunyaev-Zel'dovich effect from 
galaxy clusters, < Asz < 2 (Spergel et al. 2007). 

4.4 Additional data sets 

We also present results including constraints from the 
Union II SNIa data set of Suzuki et al. (2012), the SHOES 
program of Riess et al. (2011), and the BAO analyses of 
Percival et al. (2010), at two intermediate redshifts, and 
Reid et al. (2012), at a higher redshift. 

The SNIa data set consists of a compilation of 580 SNIa 
from a variety of sources. For the likelihood analysis of these 
data we use the COSMOMC module 16 of Suzuki et al. (2012), 
including their treatment of the systematic errors. 

For the BAO analysis, we use the results and methods 
of Percival et al. (2010), based on 2dFGRS (Colless et al. 
2003) and SDSS-II DR7 (Abazajian et al. 2009) data, by 
including bivariate Gaussian constraints on the ratio d z = 
r a (zd) I D v (z) at z = 0.2 and z — 0.35, with the corre- 
sponding covariance between do. 2 and do. 35, where r s (zd) 
is the sound horizon at the drag epoch 17 and D v (z) = 

14 http://lambda.gsfc.nasa.gov/ 

15 Here we use the five-year WMAP data (Dunkley et al. 2009; 
Komatsu et al. 2009, and companion papers). For various of our 
results we have incorporated the galaxy data by importance sam- 
pling MCMC chains that had clusters and WMAP5 data. Note 
that, for ACDM and miCDM models, the constraints on f2 m /t 2 , 
f2 c /i 2 (see comments on this parameter in Section 5.2), w, Ho and 
erg from WMAP7 data alone do not differ significantly from those 
of WMAP5 (a maximum of 15 per cent in the errors of these 
parameters and much less in most cases), and even less when 
combined with additional data sets as those in Section 4.4. We 
therefore expect a relatively small impact on the results from us- 
ing WMAP7 instead. The new WMAP9 data, which appeared in 
the last stage of the present work, have up to 36 per cent better 
errors than WMAP5 for those parameters, and promises then a 
somewhat larger impact on the results. 

16 http:/ /supernova. lbl.gov/Union/ 

17 To calculate we use the exact expression (see e.g. ap- 
pendix B of Hamann et al. 2010) instead of the approximate fit- 
ting formula of Eisenstein & Hu (1998). Note, however, that since 
Percival et al. (2010) and Reid et al. (2012) used the latter for- 
mula to fit the data, our results of the sound horizon need to be 
appropriately rescaled (see again Hamann et al. 2010). Besides 
being more accurate, the exact calculation of is independent 
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Figure 1. 68.3 and 95.4 per cent confidence contours in the erg, 7 plane for the 7+ACDM model. In the left panel, results are shown for 
the following individual data sets: galaxy growth (gal; green contours), CMB (blue contours) and cluster growth (cl; red contours). In the 
right panel, the same results are shown (for comparison purposes) together with the following combinations: cl+CMB (purple contours), 
cl+gal (magenta contours) and CMB+gal (turquoise contours). The combination of the three data sets breaks the degeneracies between 
7 and as and provides tight constraints on this plane (inner, gold contours). In the right panel, results are also shown for the combination 
gal+BAO (pale green contours). The horizontal, dashed lines mark 7 = 0.55 (GR). 



[(1 + z) 2 D A (z) 2 cz/H(z)] 1/a . Percival et al. (2010) showed 
that these results can also be recast as approximately in- 
dependent Gaussian constraints on do. 275 and the ratio of 
the distance scales D v (0.35)/D v (0.2). For the BAO results 
of Reid et al. (2012), based on CMASS BOSS data, we ex- 
tend the bivariate Gaussian likelihood of Section 4.2 to a 
trivariate Gaussian likelihood by including a constraint on 
a = [d z ]fiduciai/d z = 1.023 ± 0.019, at z = 0.57, and the 
corresponding correlation coefficients r/o- sa = —0.0086 and 
r Fa = -0.080 (see Section 6.4 of Reid et al. 2012). 

Note that the overlap between the ranges in redshift of 
the SDSS-II DR7 luminous red galaxy sample (0.16 < z < 
0.47) used in the Percival et al. (2010) results and the SDSS- 
III DR9 CMASS sample (0.43 < z < 0.7) used in the BOSS 
results (Reid et al. 2012) is very small (see this comparison 
e.g. in Anderson et al. 2012). Therefore, we assume that the 
two BAO measurements are essentially independent and can 
be straightforwardly combined. 

For the SHOES analysis, we use a Gaussian prior on 
Ho = 73.8 ± 2.4 km s Mpc . This measurement is based 
on Hubble Space Telescope optical and infrared data for over 
600 Cepheid variables in the host galaxies of 8 nearby SNIa 
(Riess et al. 2011). 



of the standard assumptions used to obtain the fitting formula, 
and therefore valid for other models. Interestingly, though, for 
the extended models used here we find no significant differences 
in the results obtained from using either calculation. 



4.5 MCMC implementation 

To calculate the parameter posterior probability distribution 
functions (pdf's) we use the Metropolis-Hastings Markov 
Chain Monte Carlo (MCMC) algorithm, as implemented in 
the code CosmoMC 18 (Lewis & Bridle 2002). We employ a 
modified version of this code that includes additional mod- 
ules for the likelihood analyses of the cluster growth exper- 
iment (MlOa) and the / gaB experiment (Rapetti et al. 2005; 
Allen et al. 2008) 19 . In this version of the code, we have also 
incorporated the RSD, AP effect and BAO (BOSS) analyses 
as CosmoMC modules 20 . We also use a modified version of 
the code CAMB 21 (Lewis et al. 2000) that includes 7 in the 
analysis of the ISW effect of the CMB data (Rapetti et al. 
2009). 

For our most general model, we simultaneously fit a 
total of 34 parameters. From these, 8 are cosmological pa- 
rameters and 26 are used to model astrophysical variables 
and marginalize over systematic uncertainties: 1 for CMB 
(see Section 4.3), 7 for / gas (see details in Allen et al. 2008) 
and 18 for cluster growth data (see Sections 2.2, 4.1.1, and 
M10a,b for full details). 



18 http:/ /cosmologist.info/cosmomc/ 

19 http:/ /www. slac.stanford.edu/~drapctti/fgas_module/ 

20 We use these and standard modules for the other data sets in 
Section 4.4 cither to run new MCMC chains or to perform impor- 
tance sampling on existing ones. For selected examples, we have 
explicitly checked that both methods provide the same results. 

21 http://camb.info/ 
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Table 1. Marginalized mean values and 68.3 per cent confidence limits for the 7+ACDM and 7+uiCDM models using various subsets 
of the data. When the CMB data is not included the BBNS and SHOES priors arc used. For gal+BAO and 7+ACDM, a slightly tighter 
constraint than the prior on Hq is obtained. 
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For analyses without CMB data, we fix n s to 0.95 
since, for such analyses, n s is degenerate with erg (see 
MlOa). For these analyses, we also use Gaussian priors 
on Ho from the SHOES program (Riess et al. 2011), and 
n b h 2 = 0.0213 ± 0.0010 from Big Bang Nucleosynthesis 
(BBNS) studies (Pettini et al. 2008). 



5 RESULTS 

5.1 Constraints on the 7+ACDM model 

The left panel of Figure 1 shows the joint constraints in the 
crg,7 plane for the 7+ACDM model. The green contours 
show the constraints obtained from the RSD and AP effect 
data from WiggleZ, 6dFGS and CMASS BOSS (hereafter 
referred as galaxy/gal data); the blue contours those from 
the CMB data; and the red contours those from the cluster 
abundance and /g aB data (hereafter referred as cluster/cl 
data). Combining the cluster+CMB+galaxy data we obtain 
the tight constraints shown by the gold contours. 

As shown in the figure, individually, the CMB and 



galaxy data exhibit significant degeneracies in the 08,7 
plane, as expected (see Section 3). For the cluster data, the 
correlation between these two parameters is much weaker, 
enabling independent constraints on both parameters. 22 Im- 
portantly, the constraints from the three independent exper- 
iments (which are affected by very different systematic un- 
certainties) are in excellent agreement. This agreement mo- 
tivates us to combine the constraints, leading to the results 
shown in the inner, gold contours. Combining the three data 
sets we obtain marginalized constraints on 7 = 0.5701q'°63 
(in good agreement with GR) and as — 0.785 ± 0.019 (see 
also Table l). 23 



22 When combining gal+BAO we also obtain results that are 
more comparable to those of the clusters due to the degeneracies 
broken by this combination (sec further details in the text). 

23 For this paper, we quote marginalized mean values and cen- 
tral credible intervals. For the latter, equal fractions of the vol- 
ume of the posterior lie on each side of the interval (see e.g. 
Harnann et al. 2007). Instead, in R10 results were presented us- 
ing marginalized peak values and minimal credible intervals, for 
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Figure 2. 68.3 and 95.4 per cent confidence contours in the w, 7 plane for the 7+toCDM model from cl+CMB+gal (gold contours 
in both panels) and the following combinations of data (right panel): cl+CMB (purple contours), cl+gal (magenta contours), and 
CMB+gal (turquoise contours). The platinum contours in the left panel correspond to adding SNIa+SHOES+BAO to the combination 
of cl+CMB+gal. The horizontal, dashed lines mark 7 = 0.55, the growth history for GR. The vertical, dot-dashed lines mark w = — 1, 
the expansion history for ACDM. This figure shows that the results are simultaneously consistent with GR and ACDM. 



If we also include SNIa, BAO and the SHOES measure- 
ment of Ho, the constraints on the growth parameters are, 
as expected, almost the same (see Table 1) although inter- 
estingly we obtain a small 4 per cent improvement in the er- 
ror in 7. For this combination we also obtain improved, tight 
constraints on the expansion parameters f2 m = 0.277 + 0.011 
and Ho = 70.2 ± 1.0 km s _1 Mpc _1 . It is worth noting that 
the addition of the BAO data alone provides almost the 
same improved constraints on the r2 m ,i/o plane as those 
from adding all three data sets (see Table 1). 

The right panel of Figure 1 shows a zoom into the cen- 
tral regions of the constraints shown in the left panel, to- 
gether with the constraints for the combinations of cl+CMB 
data (purple contours), cl+gal data (magenta contours) and 
CMB+gal data (turquoise contours). The gold, tightest con- 
tours correspond again to the combination of the three 
data sets, cl+CMB+gal. Notably, the nearly orthogonal de- 
generacies of the CMB (blue contours) and galaxy (green 
contours) constraints allow their combination (turquoise 
contours) to provide tight marginalized constraints in the 
growth plane. The area enclosed by the 95.4 per cent confi- 
dence contour in the a"8,7 plane is only slightly more than 
one third larger for CMB+gal than for the three data sets 
combined. 



which the size of the interval is minimized. For approximately 
symmetric posteriors, such as those for our combined data, both 
choices provide similar results, although by construction those 
from the former tend to be slightly more conservative. For our 
individual data sets, for which the posteriors are less symmetric, 
we show full marginalized distributions in Figures 4 and 5. 



As found by R10 24 , for the cl+CMB data (purple con- 
tours) erg and 7 are highly correlated, with a correlation co- 
efficient p = —0.85. The addition of the galaxy data breaks 
this degeneracy. With respect to constraints obtained from 
cl+CMB, those for the cl+CMB+gal provide more than a 
factor 2 reduction in the area enclosed by the 95.4 per cent 
confidence contour in the o"8,7 plane. 

In the right panel of the figure, we also show the con- 
straints from the combination gal+BAO (pale green con- 
tours), for which both data sets come from the analysis 
of different properties of galaxy redshift surveys. Interest- 
ingly, even though the baryon acoustic oscillation data on 
their own provide only constraints on expansion parameters, 
those on f2 m help in reducing the large degeneracies that 
the galaxy growth data has in the fi m , 7 and f2 m , as planes, 
with correlation coefficients of p = 0.83 and p = 0.74 for 
each plane. Adding BAO to gal we obtain then a significant 
improvement in the constraints on the growth plane o"8,7- 

5.2 Constraints on the 7+iuCDM model 

The left panel of Figure 2 shows the joint constraints in 
the w, 7 plane for the 7+wCDM model. For the combina- 
tion of our primary data sets, cl+CMB+gal, we obtain the 
gold contours. For these, we find marginalized constraints 
on w = -0.9501^ and 7 = 0.533 ± 0.080 at the 68.3 per 



Note that the results presented in R10 were for a combination 
of cluster+CMB+SNIa+BAO data. However, the constraints on 
the as, 7 plane were primarily driven by the cluster+CMB data. 
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Figure 3. 68.3 and 95.4 per cent confidence contours in the growth plane erg, 7 (left panel) and the expansion planes f2 m ,ui (middle 
panel) and Cl m , Hq (right panel) for the 7+«>CDM model. The three panels show constraints from the following combinations of data sets: 
cl+CMB (purple contours), cl+gal (magenta contours), CMB+gal (turquoise contours), and cl+CMB+gal (gold contours). The platinum 
contours correspond to adding SNIa+SHOES+BAO to the latter combination of data. The horizontal, dashed lines mark 7 = 0.55 (GR; 
left panel) and w = —1 (ACDM; middle panel). 



cent confidence level. These results are simultaneously con- 
sistent with GR and ACDM. The platinum contours in this 
panel show the joint constraints on the w, 7 plane when we 
add SNIa+SHOES+BAO to the cl+CMB+gal data 25 . In this 
case, we find marginalized constraints of w — — 0.9871q'q53 
and 7 = 0.604 + 0.078. Again the results are consistent with 
GR+ACDM. 

In the right panel of the figure, the purple contours 
correspond to cl+CMB, the magenta contours to cl+gal, the 
turquoise contours to CMB+gal and the gold contours again 
to the combination of the three data sets. The horizontal, 
dashed and vertical, dot-dashed lines mark 7 = 0.55 (GR) 
and w = — 1 (ACDM), respectively. 

Comparing the cl+CMB with the cl+CMB+gal results, 
we find 46 and 62 per cent improvements in the constraints 
on 7 and erg. It is also worth noting that the improvement in 
the joint measurement of w and 7 is larger than that for each 
individual parameter. We find more than a factor 3 reduc- 
tion in the area enclosed by the 95.4 per cent confidence con- 
tour of the joint w, 7 constraints. Note that the correlation 
between w and 7 increases from p = —0.47, for cl+CMB, to 
p = —0.66, for cl+CMB+gal, which suggests that additional 
constraints on w might also help improving those on 7. In 
fact, even though SNIa and SHOES data provide direct addi- 
tional constraints on only cosmic expansion parameters, for 
which we obtain e.g. a 27 per cent improvement on w when 
adding them to cl+CMB+gal, the combined, marginalized 
constraints on 7 represent a small improvement of 4 per cent 
due to the correlation between w and 7. For these data sets 
combined, cl+CMB+gal+SNIa+SHOES, the correlation in 
the w, 7 plane is still of p = —0.65. Interestingly, the corre- 
lation between 7 and ft c h 2 is also relatively large, p = 0.72, 26 

25 From Table 1, note that the main additional constraint on 
this plane comes from the SNIa data. The SHOES and BAO data, 
though, significantly help in constraining the combination of pa- 
rameters f2 m and Hq. 

26 The correlation between 7 and the CMB shift parameter, 



which indicates that e.g. the significant improvements in the 
constraints on this parameter from the CMB measurements 
of the Planck satellite 27 should help with constraining 7. 

Figure 3 shows constraints for the same model and sub- 
sets of the data for three different planes: the growth plane 
erg, 7 (left panel) and the expansion planes Q m , w (middle 
panel) and O m , Ho (right panel) . The left panel of this fig- 
ure shows that the correlation between erg and 7 reduces 
dramatically from cl+CMB (purple contours), p — —0.56, 
to cl+CMB+gal (gold contours), a negligible p = 0.08. The 
reduction in the area enclosed by the 95.4 confidence con- 
tour in this growth plane when adding gal to the cl+CMB 
data is substantial. 

The platinum contours in the three panels correspond to 
the constraints obtained when adding SNIa+SHOES+BAO 
to the cl+CMB+gal data. The improvement in the growth 
plane of the left panel of the figure is small while those in 
the expansion planes of the middle and right panels are sig- 
nificant due to the degeneracy breaking power of the addi- 
tional data in these planes. For this model, the combined 
constraints on O m = O.278t[5:oii an d H = 70.0 ± 1.3 km 
s _1 Mpc _1 are again very tight. 



6 DISCUSSION 

6.1 Comparing results 

For a ACDM expansion model, and combining galaxy and 
CMB data, recent studies have presented constraints on 7 

which is highly correlated with Q, c h? (p = 0.92) and is purely 
an expansion parameter at high- z (see the definition and further 
details in e.g. Komatsu et al. 2009), is also large, p = 0.68. The 
correlations between 7 and f2 m (= C^+^c), p = 0.32, 7 and Q^h 2 , 
p = 0.11, and 7 and Ho, P = 0.13, are significantly smaller. As 
shown in the right panel of Figure 3 for other data combinations, 
the correlation between Q m and Ho is also large, p = —0.82. 
27 http://www.esa.int/Planck 
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Figure 4. Marginalized probability distribution functions for 7 (left panel) and ag (right panel) for the 7+ACDM model. Results are 
shown for the following data sets: CMB (blue, dotted line), gal (green, long-dashed line), el (red, dashed line) and cl+CMB+gal (gold, 
solid line). In the left panel, the vertical, dashed line marks 7 = 0.55 (GR). 



that are similar to and in agreement with ours. For exam- 
ple, Hudson & Turnbull (2012) combined data from two pe- 
culiar velocity surveys at low redshifts (Davis et al. 2011; 
Turnbull et al. 2012) and RSD (but not AP effect) data from 
various galaxy surveys. Samushia et al. (2013) used primar- 
ily RSD, AP effect, and BAO data from the CMASS BOSS 
results of Reid et al. (2012) together with RSD and AP ef- 
fect data from other surveys. In their combined results, both 
studies include WiggleZ, 6dFGS and CMASS BOSS data, as 
we do here, in addition to other galaxy and expansion data 
sets. Both analyses use CMB data from WMAFJ. The for- 
mer study uses previous results from CMB and expansion 
data only as a prior, while the latter uses the full CMB 
likelihood 28 . Neither of these analyses, however, use the low 
multipoles of the CMB to constrain 7 with the ISW effect 
(see Section 6.1.1). 

Note also that these studies include BAO constraints 
from Percival et al. (2010) and Reid et al. (2012), respec- 
tively. As discussed in Section 6.1.2, both BAO data sets 
(and especially the latter) prefer larger values for fl m , which 
in combination with growth data implies a preference for 
larger values of 7. This, together with the fact that these 
works do not include the cluster data or the ISW effect 
constraints from the CMB data, which both prefer smaller 
values for 7, is consistent with their results on 7 being at 
the high end of ours in Table 1. Although all these results 
and those in Table 1 are consistent with GR (7 ~ 0.55), 
the differences highlight the importance of studying each 



28 Note that it is important to include 7 in the full CMB analysis, 
in combination with the other experiments, to account for all the 
degeneracies of the CMB parameters with both expansion and 
growth parameters, such as e.g. that of 7 with f! c h 2 or the CMB 
shift parameter. If these covariances are not included, one may 
obtain spuriously tight results. 



individual data set as well as their various combinations in 
detail before combining all of them. For upcoming, more 
statistically powerful data sets, this will also be increasingly 
important. 

6.1.1 ISW effect 

Even though the ISW effect has only a relatively small im- 
pact on the combined results, it is not negligible. Using 
our analysis, it is interesting to compare results including 
or not the ISW effect for 7. For the 7+ACDM model and 
the combination CMB+gal, we obtain 7 = 0.607^°' oso 

with- 
out the ISW effect, which as expected (see Section 6.2) is 
slightly higher (3 per cent) than our default result (see Ta- 
ble 1) and weaker by 8 per cent. 29 For cl+CMB, we obtain 
7 = 0.432±g:i|§ and a s = 0.842 ± 0.057, which are 20 and 
16 per cent weaker than our default results (see also a sim- 
ilar comparison in Rapetti et al. 2009). For cl+CMB+gal, 
7 = 0.585 ± 0.067 is only 6 per cent weaker than the corre- 
sponding result including the ISW effect for 7. 

6.1.2 Adding the BAO data 

Our results show (see both Table 1 and the left panel of Fig- 
ure 4) that, compared with the cluster and CMB data, the 

29 For the same data but for the 7+uiCDM model, we obtain 
7 = 0.547^Q gg 2 , higher by 4 per cent and weaker by 5 per cent, 
and w = —0.914+0.073, tighter by 10 per cent. For this model, the 
background expansion parameter ui modifies the ISW effect in an 
approximately opposite way to that of the density perturbation 
parameter 7. Therefore, by artificially ignoring 7 in the calcu- 
lation of the ISW effect, tighter constraints on w are obtained. 
This is comparable to what happens for a ui-model of a dark 
energy fluid when its dark energy perturbations are erroneously 
not taken into account (Weller &: Lewis 2003; Bean & Dore 2004; 
Rapetti et al. 2005). 
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Figure 5. Marginalized probability distribution functions for 7 (left panel) and w (right panel) for the 7+uiCDM model. Results 
are shown for the following data sets: cl (red, dashed line), cl+CMB (purple, dotted line), cl+CMB+gal (gold, solid line) and 
cl+CMB+gal+SNIa+SH0ES+BAO (platinum, solid-thin line). The vertical, dashed lines mark 7 = 0.55 (GR; left panel) and w = -1 
(ACDM; right panel). Note that 7 (cosmic growth) and w (cosmic expansion) are measured simultaneously with similar precision. 



combination gal+BAO prefers significantly larger values of 
f2 m , and therefore of 7 due to the covariances between fi m 
and 7, and f2 m and as (see Section 5.1). 30 Also, for any of 
the data set combinations in Table 1, the addition of the 
BAO data shifts the constraints on fi m and 7 to larger val- 
ues. Adding BAO to all the other data sets combined and for 
the 7+U1CDM model, we have increases of 9 and 12 per cent 
for each parameter 31 . It is interesting to note, though, that 
using only the BOSS BAO data set, we obtain similar shifts 
of 7 and 11 per cent, although slightly weaker constraints 
on O m , and similar constraints on 7. Using instead only the 
BAO data set of Percival et al. (2010), we find about half of 
those increases, 5 per cent for both parameters, and also a 
bit weaker constraints on fi m . The constraints on 7, though, 
are slightly tighter due to the reduction in the tension with 
the other data sets. The mild tension on ri m between the 
BAO and the other data sets translates in some cases into a 
smaller constraining power for 7 (and also for w) when com- 
bining them. Table 1 shows e.g. that for 7+uiCDM, adding 
both BAO data sets to cl+CMB+gal provides somewhat 
weaker constraints on 7, and also that these are 13 per cent 
weaker than those for instead adding SNIa to cl+CMB+gal. 
In addition, using all the data sets combined except BAO, 
we obtain the tightest constraints on 7 for 7+ACDM, and 
on both 7 and w for 7+wCDM. However, the increase in 
constraining power on these parameters is small compared 



30 For 7+ACDM, using the CMB data alone we have f2 m = 
0.260 ± 0.030. For gal+BAO, using only the BAO data set from 
Percival et al. (2010), we obtain U m = 0.345 ± 0.050 and 7 = 
0.719ig24g, and using instead only the BOSS BAO data set, 
C m = 0.417±0.073 and 7 = 0.978±q 3 g° , which are clearly larger. 

31 Note that the shift between the gold, solid and platinum, solid- 
thin lines in Figure 5 is mainly due to the addition of BAO. 



with the decrease in constraining power on the other expan- 
sion parameters when not using BAO. 

The BAO and SHOES data also present a mild tension in 
the direction of the well-known degeneracy between Ho and 
O m (sec e.g. Hinshaw et al. 2012). The addition of SHOES 
to any of the data combinations in Table 1 that include our 
primary data sets, shifts Ho to larger values, and therefore 
f2 m to smaller values through the correlation between these 
two parameters. 



6.2 Constraining power 

As discussed in Section 5.2, the combination CMB+gal pro- 
vides tight constraints on the <7s,7 plane (see Figure 1) due 
to the complementarity between the constraints from the 
individual data sets. However, the large degeneracies of the 
individual constraints make the combination prone to po- 
tential biases from systematic uncertainties. The left panel 
of Figure 4 shows that for galaxies alone (green, long-dashed 
line) the marginalized pdf for 7 has a large tail toward values 
higher than that for GR, although interestingly the peak is 
close to the GR value (vertical, dashed line). 32 On the other 
hand, for the CMB (blue, dotted line) values larger than 
that for GR are significantly constrained by the data (due 
to the ISW effect), while lower values are largely uncon- 
strained and degenerate with as, which has an extended tail 
toward large values (see the right panel of the figure). 

From comparing the normalized pdf 's in the figure, it is 
worth noting that while the constraining power of the clus- 
ter data on 7 (left panel) and as (right panel) is notably 



32 Adding the expansion data set BAO to gal shortens this tail 
(see Figure 1) and shifts the peak to a larger value (see Table 1). 
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better than that of the CMB or galaxy data, the combina- 
tion cl+CMB+gal is much more powerful than the cluster 
data alone. Note also that the power of the current data for 
constraining erg (right panel) is considerably greater than for 
constraining 7 (left panel). 

6. 2. 1 Full model: 7 +w CDM 

For our most general model, 7+wCDM, only the cluster 
data can alone constrain this model at a significant level. 
We obtain w = -1.021+2;^ and 7 = 0.507t%$a ( see also 
Figure 5). 33 Since our other primary data sets do not have 
strong direct constraints on 7 (see Section 3), their con- 
straining power depends critically on the complexity of the 
model used. For our extended model, we allow departures 
from the standard expansion and growth histories equally. 
Combining all our data sets, we obtain the tightest and most 
robust results to date on this model. The addition of SNIa, 
SHOES and BAO data is particularly helpful for constraining 
the expansion parameters in this model. The right panel of 
Figure 5 shows that when we include these data sets (plat- 
inum, solid-thin line) the constraining power on w clearly 
increases. The figure also shows the progression in the pdf 's 
of 7 (left panel) and w (right panel) when adding one at a 
time the other primary data sets to the cluster data. Re- 
markably, for these combinations (as well as for the others 
of the primary data sets) we can measure at the same time 
7 (cosmic growth) and w (cosmic expansion) with similar 
precision. 



7 CONCLUSIONS 

We have combined cluster growth and / gas data from RASS 
and CXO, CMB data from WMAP, and RSD and AP ef- 
fect data from WiggleZ, 6dFGS and CMASS BOSS to si- 
multaneously constrain the evolution of cosmic structure 
and background expansion. To test for consistency with GR 
and ACDM, we have used convenient parameterizations: 
f2 m , Ho and w for the expansion history, and erg and 7 for 
the growth history. We find that the combination of clus- 
ters+CMB+galaxies breaks key degeneracies in the growth 
plane, erg versus 7, for the data sets individually. In combi- 
nation, the data provide tight, robust constraints that are 
in excellent agreement with GR+ACDM. 

Fixing w — —1, we obtain marginalized constraints 
on the growth parameters erg = 0.785 ± 0.019 and 7 = 
°- 570 ±ao63- Including SNIa, SHOES and BAO data we ob- 
tain 7 = 0.616 ± 0.061. Allowing w to vary, we have erg = 
0.780 ± 0.020 and 7 = 0.533 ± 0.080 for the combination of 
clusters+CMB+galaxies. For this, we find a correlation be- 
tween w and 7 of p = -0.66. Including SNIa+SH0ES+BAO, 
we obtain fi m = 0.278±°;on> H = 70.0 ± 1.3 km s^Mpc" 1 
and w = —0.987^0 053 f° r the expansion parameters, and 
erg = 0.789 ± 0.019 and 7 = 0.604 ± 0.078 for the growth 
parameters. 



33 For the combination gal+BAO (sec Tabic 1), we obtain similar 
constraints on w and 7, while those on fi m and erg are notably 
weaker. 



Our results highlight the potential of combining forth- 
coming galaxy cluster data (from e.g. the South Pole Tele- 
scope [SPT], the Atacama Cosmology Telescope [ACT], 
yLMM-Newton wide area surveys, the Dark Energy Survey 
[DES], the Large Synoptic Survey Telescope [LSST], and 
the extended ROentgen Survey with an Imaging Telescope 
Array [eROSITA]), CMB data (from e.g. SPT, ACT and 
Planck), and galaxy data (from e.g. SDSS-III, the Subaru 
Measurement of Images and Redshifts [SuMIRe] project, 
BigBOSS, the Dark Energy Spectrometer [DESpec], and the 
Euclid mission) for constraining dark energy and modified 
gravity models. 
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